Computational study of (111) epitaxially strained ferroelectric perovskites BaTiOa 

and PbTi0 3 
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The phase transition behaviour of PbTi03 and BaTi03 under (111) epitaxial strain is investigated 
using density-functional theory calculations. Prom tensile strains of +0.015 to compressive strains of 
-0.015, PbTiOs undergoes phase transitions from C2 through two Cm phases and then to R3m. The 
total polarisation is found to be almost independent of strain. For the same range of strains BaTi03 
undergoes phase transitions from a single Cm phase, through R3m and then to R3m. In this case 
the application of compressive strain inhibits and then completely suppresses the polarisation on 
transition to the non-polar R3m phase. 

PACS numbers: 77.55. +f,77.84.Dy,77.22.Ej,77.80.Bh 



I. INTRODUCTION 

Perovskite materials have been extensively studied due 
to their numerous structures and properties^. At room 
temperature bulk BaTiC>3 and PbTiC>3 are ferroelectric 
with spontaneous polarizations of around 25 ^Ccm~ 2 and 
80 /iCcm -2 , respectively. Epitaxial strain can have a 
considerable effect on the structural and electrical prop- 
erties of perovskite thin films 2 -^^. The properties of 
thin films can, therefore, be fine-tuned by growing the 
films epitaxially and coherently on a substrate with a 
desired lattice constant. 

The effect of strain has already been studied for per- 
ovskite thin films oriented in the (001) direction on a 
cubic substrate^^^. (Ill) epitaxial films have differ- 
ent symmetries than (001) films and, therefore, the be- 
haviour will also be different. The temperature-strain 
phase diagram of PbTiC>3 (111) films has been modelled 
using a thermodynamic mean-field approach^. For room 
temperature and below, it was predicted that the film 
should have 3m point symmetry for compressive strains 
and two different m symmetry phases for tensile strains, 
with first-order transitions between the phases. In this 
paper we investigate the stability of various phases of 
(111) films of PbTiC>3 and BaTiC>3 under epitaxial strain 
at T=0 K using density-functional theory calculations. 
We show that the zero-temperature phases of PbTiC>3 
mostly agree with the higher-temperature results of Ref. 
HcL BaTiOs exhibits a qualitatively different behaviour 
from PbTiC>3 and, surprisingly, the application of com- 
pressive strain is found to suppress the polarization. 



II. METHOD 



rameters to be (1,0,0) and (— 5, 2 ,0) and scaling both 
lattice vectors with mismatch strain. We have not in- 
cluded a surface as a ferroelectric field in vacuum gives 
rise to a depolarizing field, which suppresses polarization. 
This approach allows us to separate the effect of epitaxial 
strain from surface effects. 

The ISOTROPY code^i and the Bilbao Crystallo- 
graphic Server— were used to find space groups com- 
patible with a perovskite (111) film. Fig. Q] shows the 
various subgroups of R3m with a possible 5-atom prim- 
itive cell and their symmetry relations. The monoclinic 
geometry of the substrate implies that PI will actually 
have a higher C2/m symmetry, and R3 actually has the 
higher R3m symmetry for perovskites as the mirror-plane 
breaking Wyckoff positions are not occupied. Therefore, 
eight possible space groups exist and are shown in Table[H 
The possible unstrained bulk symmetries that reduce to 
these space groups when (111) strain is applied are listed 
in Table Hlj with the Cm phase divided in three types as 
shown in Fig. [21 

R-3m 



R3m 




FIG. 1: (Color online) The subgroups of R3m with a possible 
5-atom primitive cell. Crossed symmetries are not present in 
(111) epitaxial perovskite. 



The parent space group for a (111) epitaxially strained 
perovskite film is Rim. In this work we used 15-atom 
unit cells, but only considered symmetries compatible 
with a 5-atom primitive cell. Periodic boundary con- 
ditions were used, and the effect of the substrate is im- 
plicitly included by constraining the a and b lattice pa- 



For the (111) perovskite, fewer polarization directions 
are allowed than in the (001) case. The non-polar R32 
and C2/m space groups are not expected to be seen in the 
(111) perovskite^, and they have not been considered. 
For displacements strictly perpendicular to the plane, the 
polarization P3Z and the space group Rim is obtained. 
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TABLE I: Description of the non-polar and polar (111) epi- 
taxial perovskite symmetry groups with a possible 5-atom 
primitive cell. Vectors are given in the standard rectangular 
Cartesian basis (a, 2 t^, %)• Np is the number of free param- 
eters in atomic fractional coordinates when we restrict the 
primitive cell to 5 atoms, although the symmetries marked 
with * would allow a 15-atom primitive cell, with more free 
parameters in fractional coordinates. 
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TABLE II: Considered strained symmetries and bulk symme- 
tries that reduce to them when (111) strain is applied. 



a + 2b 




FIG. 2: Possible polarization directions in the mirror plane of 
the Cm symmetry. Tilde denotes the presence of symmetry 
at zero strain. 



R3m is predicted to be stable for compressive strains in 
PbTiOy^ -. i?3 is a subgroup of Rim that has polariza- 
tion in the same direction but additional displacements 



of atoms to break the mirror plane. 

For displacements parallel to a lattice vector in the 
plane, the polarization P\a and the space group C2 is 
obtained. The stability of such a PbTiC-3 phase is consid- 
ered uncertain in Ref. [HI as it depends on the selection 
of model parameters. 

For displacements parallel to a + 2b (or, equivalently, 
a — b or a + 2b) in the plane, the space group Cm is 
obtained, in which the polarization is only constrained to 
lie in a mirror plane perpendicular to the epitaxial plane. 
According to the Landau theory, Cm is stable for tensile 
strains in PbTiOjr^ (when using the parameters from 
Ref. [13), and the transition Cm — > Rim is first-order. 
Due to symmetry considerations, transition Cm — > C2 
also has to be first-order. 

Within the monoclinic Cm space group, there are three 
inequivalent possible energy minima, with isotropic phase 
transitions between them. They correspond to the bulk 
space groups Rim, Amm2 and PAmm and pseudocu- 
bic bulk polarization directions (111), (110) and (100), 
respectively-^. In Ref. 0, the stability of the monoclinic 
phase in bulk is studied with high-order Devonshire the- 
ory, and possible monoclinic Ma and Mb phases (along 
with Mc, which would have Pm symmetry that is not 
possible in a (111) biaxially strained perovskite) originate 
in the eighth-order expansion of the thermodynamic po- 
tential. Our Amm2 minimum corresponds to Mb and 
PAmm minimum to Ma, as displayed in Fig.[2l 

Finally, a triclinic phase with no symmetry constraints 
(PI) would require twelfth-order terms in Ref. and is 
predicted to be unlikely in the bulk. No such phase is 
considered in Ref. [l^ either, and it is left outside the 
present study. 

In a strained film, there are nonzero stresses <j\ and 
o"2 (in the standard Cartesian axes) in the epitaxial di- 
rection. Therefore, the thermodynamic potential that is 
minimized is the elastic enthalpy 



H = U + e m (a 1 + cr 2 ), 



(1) 



where U is the internal energy and e m the in-plane 
strain relative to the relaxed lattice constant of the phase 
in question. The phase with lowest elastic enthalpy at 
given mismatch strain is the equilibrium state of a co- 
herent epitaxial single-domain film. The internal energies 
and stresses for each phase as a function of strain are ob- 
tained from density-functional theory (DFT) calculations 
after relaxing atomic positions and unit cell parameters 
within the constraints of Table [H 

The DFT calculations were performed using the Vi- 
enna ab-initio simulation package (VASP) 15 ' 16 . The core 
states were represented using the PAW metho d 17 ' 18 and 
the semicore Pb, Ba, and Ti were treated as valence elec- 
trons. The planewave cutoff energy was 700 eV. The ex- 
change and correlation energy was described using the lo- 
cal density approximation (LDA). The polarization of the 
phases was calculated using Berry phase formalis m 19 ' 20 . 
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The hexagonal (monoclinic for Cm) conventional cell 
contained 15 atoms and consisted of three identical 5- 
atom primitive cells. A T-centred Brillouin zone mesh of 
6x6x4 was used. 

Ionic relaxations were performed for various c/a ratios 
(for all space groups) and c axis tilt angles (for Cm) 
to obtain the optimal cell shape and volume for each 
fixed in-plane strain. Relaxations were stopped when the 
Hellmann-Feynman forces on the atoms were below 0.01 
eV/A. 



III. RESULTS 

The relaxed DFT bulk cubic lattice values were 3.895 A 
for PbTi0 3 and 3.955 A for BaTi0 3 , and they were used 
as the unstrained lattice constants. The calculated c/a 
ratios in the bulk tetragonal phase were 1.04 for PbTiC>3 
and 1.003 for BaTiC>3, and corresponding polarizations 
in the bulk tetragonal phase were 0.77 Cm -2 for PbTi0 3 
and 0.21 Cm -2 for BaTiC>3. As expected, these values are 
somewhat smaller than the experimental values, which is 
typical of the LDA approximation. The experimental 
cubic lattice constants are 3.97 A2i and 4.00 A2 2 ., the 
tetragonal c/a ratios are 1.07 and 1.01 and the polariza- 
tions are 0.80 Cm" 2 ^ and 0.26 Cm" 224 , for PbTi0 3 and 
BaTiC>3, respectively. 

The internal energies and corresponding elastic en- 
thalpies for the considered epitaxial phases, with error es- 
timates, are displayed as a function of strain in Fig.[4]for 
BaTiC>3 and Fig. [5] for PbTiC>3. In the enthalpy curves, 
the energetically favored phases are denoted. The error 
in enthalpy is proportional to relative strain squared. For 
PbTiC>3, two separate but almost identical energy curves 
for the Cm phases were obtained, corresponding to dif- 
ferent energy minima within the same symmetry, by re- 
laxing the unit cell starting from three different starting 
configurations. For the strains considered, -R3 was found 
to relax into the R3m symmetry for both PbTiC>3 and 
BaTiC>3, so no separate R3 curve is plotted. 

In-plane (Pi and P2), perpendicular to plane (-P3), and 
total (P) polarization of the perovskites as a function of 
strain are displayed in Fig. [4] for BaTiC>3 and Fig. [5] for 
PbTiC>3. The components displayed are strictly from the 
phase of lowest enthalpy for each strain as denoted, al- 
though the error due to discretization in strain, c/a ratio 
and c axis angle makes the preferred phase uncertain at 
large compressive or tensile strains and between the two 
PbTi0 3 Cm phases. 



A. BaTiO :j 

For BaTiC>3, only a single energy curve for the Cm 
phase could be obtained, with relaxation to a single en- 
ergy minimum instead of transitions between them. This 
energy minimum most closely corresponds to a nearly 
R3m structure, with polarization towards a Ba atom as 



indicated in Fig. [3l and virtually no rotation of polar- 
ization with increasing strain. The result is in accor- 
dance with the known tendency of BaTiC>3 to adopt 
the rhombohedral phase in the bulk. The transition 
Cm — > R2>m is seen to be first-order (with polarization 
vector as the order parameter), as predicted by an ex- 
tension of Devonshire theory^. The other BaTiC>3 tran- 
sition, ferroelectric-to-paraelectric (R3m — > R3m) tran- 
sition induced by compressive strain, is similarly seen 
to be first-order. Although the polarized R3m phase 
has a lower internal energy than the R3m phase, the 
lower stresses observed in the paraelectric phase make it 
favored at compressive strains when elastic enthalpy is 
considered. 
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FIG. 3: Approximate polarization directions in the mirror 
plane of the Cm symmetry. BaTiOs adopts a single energy 
minimum: dashed arrow indicates polarization of BaTiOs 
at tensile strain 0.0025. PbTiOg adopts two energy min- 
ima: black arrows indicate polarization of PbTi03 Mb and 
Ma phases at tensile strains and 0.005, respectively. Tilde 
means symmetry direction at zero strain. 



B. PbTiOa 

In PbTiC>3, the Cm polarization direction at tensile 
strains is dominated by the tendency to adopt a nearly 
PAmm {Ma) or Amm2 (Mb) structure. These polar- 
ization directions in the mirror plane of Cm phase are 
sketched and labelled in Fig. ED No relaxation to a nearly 
R3m structure, as in BaTiC>3, is observed. Therefore, 
when the tensile strain is increased, a transition from 
Mb to Ma structure, with consequent flipping of one 
polarization component, is seen in Fig. [5j Such a strain- 
induced transition is indeed predicted for (111) PbTiC>3 
in Ref. [l(J and all the PbTiC>3 transitions are first- 
order (polarization vector changes discontinuously) as 
predicted. 

This is in apparent contradiction with the fact that 
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FIG. 4: (Color online) (Top) The internal energies (per 5- 
atom primitive cell) given by VASP for considered phases of 
BaTiOs. (Center) Elastic enthalpies (per 5-atom primitive 
cell) as given by Eq. ^ for considered phases of BaTi03. 
The error bars indicate approximate error in enthalpy due to 
stress inaccuracy in VASP and discretization in strain, c/a 
ratio and c axis angle. The vertical lines denote approxi- 
mate points of phase transition. (Bottom) Polarization of the 
BaTiC>3 phase with lowest elastic enthalpy at each strain. P is 
the total polarization and Pi, P2 and P3 are the polarization 
components in the a, and z directions, respectively. 



FIG. 5: (Color online) (Top) The internal energies (per 5- 
atom primitive cell) given by VASP for considered phases of 
PbTiOs. (Center) Elastic enthalpies (per 5-atom primitive 
cell) as given by Eq. JTJ) for considered phases of PbTiOs. 
The error bars indicate approximate error in enthalpy due to 
stress inaccuracy in VASP and discretization in strain, c/a 
ratio and c axis angle. The vertical lines denote approxi- 
mate points of phase transition. (Bottom) Polarization of the 
PbTiC>3 phase with lowest elastic enthalpy at each strain. P is 
the total polarization and Pi, P2 and P3 are the polarization 
components in the a, a+ 2- b and z directions, respectively. 



bulk PbTiC>3 is known to adopt PAmm symmetry when 
unstrained; however, the unstrained (111) epitaxial case 
is different from the bulk. Even at zero strain, the (111) 
substrate inhibits the tetragonal elongation of unit cell 
required for the PAmm phase to be stable, and the 
Amml (Mb) geometry is adopted instead. When in- 
creasing tensile strain, an Ma structure closer to PAmm 
is adopted, as it allows the polarization to align more 
parallel to the strain. 

The single enthalpy differences between the PbTiC>3 
Cm structures are smaller than possible errors in en- 



thalpy, but the enthalpy curves of Fig. [5] as a whole dis- 
play slight differences between the Ma and Mb struc- 
tures. Polarization in the two phases rotates very little 
in the narrow strain range where each phase is stable, 
though it would rotate significantly more if Cm were 
stable at larger tensile strains. At large tensile strains, 
although the Cm phase still has the lowest internal en- 
ergy, the C2 phase is favored in terms of elastic enthalpy, 
as its relaxed lattice constant is higher than that of the 
Cm phase. Therefore, large tensile strain (around 1 %) 
causes PbTiC>3 polarization to align along a lattice vector 
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TABLE III: Calculated LDA band gaps for lowest-enthalpy 
phases at different in-plane strains of (111) epitaxial per- 
ovskite. In the pseudocubic Brillouin zone, the band gaps 
are indirect T-X and F-R gaps for PbTiC>3 and BaTi03, re- 
spectively. Energies are in eV. 

in the (111) plane; the structure is similar to the Amm2 
phase in bulk, but with lower symmetry due to strain. 

C. Bandgaps 

Finally, the LDA bandgaps for the favored phases 
were calculated, to compare to LDA bandgaps of cubic 
BaTiC-3 and PbTi0 3 . While LDA seriously underesti- 
mates bandgap, any strain-induced changes in the LDA 
bandgaps may indicate widening or closing of the band 
gap due to strain. The results for the two perovskites are 
shown in Table IIIII The band structures and gaps are 
very similar to calculations on cubic perovskitesA, and 
the T-X type of indirect gap in PbTiC>3 agrees with opti- 
cal measurements of the room temperature bulk 2 ^ It can 
be seen that all non-pseudocubic phases have the effect of 
widening the gap compared to the calculated bulk cubic 
values, but there are no clear trends or major changes in 
the size of the bandgap. 

IV. DISCUSSION AND CONCLUSIONS 

BaTiC>3 has Cm symmetry for large tensile strains and 
undergoes a phase transition to the R3m phase at a strain 
of ~0.002. In both R3m and Cm phases, the total po- 
larization decreases approximately linearly as compres- 
sive strain is increased. The most surprising result is 
the adoption of the non-polar R3m phase below a strain 
of —0.004. Therefore, compressive strain appears to in- 
hibit and, finally, suppress, ferroelectricity in (111) epi- 
taxial BaTiC>3. There are observations of such inhibition 
of remanent polarization in (111) epitaxial BaTiC>3 on 
LaNiC>3-coated SrTiC>3 substrate^ 7 , which would inflict 
compressive strain due to lattice mismatch, and in (111) 
BaTiOs/SrTiOs superlattices^. However, it is also pos- 
sible that a cell-doubling transition may occur, as only 
phases with a 5-atom primitive cell were considered in 
this study. 

Both perovskites have a strain region in which the 
R3m phase is stable and undergoes a first-order phase 
transition into the Cm phase when tensile strain is in- 
creased. The geometry of the Cm phase differs in the 
two perovskites, with a single minimum in BaTiC-3 and 
first-order isomorphic transition between two Cm phases 
in PbTiC>3, with switching of the vertical polarization 



component. Such a transition was predicted in Ref. [l(] 
although their phase diagram has a lower limit of —50° C 
and our first-principles calculations are valid at absolute 
zero. Another indication of this isomorphic transition is 
the observation of inversion of the vertical polarization 
in (111) PZT thin-film capacitors^. However, the strain 
in these thin-film capacitors is unknown. 

In PbTiC>3 it is observed that as tensile strain grows, 
the C2 phase (with in-plane polarization) is adopted. In- 
deed, in Ref. [l(J it is noted that such a phase could be 
stable instead of Cm, depending on the parameters used 
in the thermodynamic expansion. The total polarization 
of the stable phases of PbTiC>3 is observed to be almost 
independent of the applied strain. The direction of polar- 
ization, but not its magnitude, changes with (111) strain. 
Such a suppressed dependence of polarization on strain 
has been observed earlier in tetragonal (100) epitaxial 
PbZro.2Tio.sO3 (PZT)2£. This was proposed to be due 
to the fact that the polarisation of PZT is already large 
so the Pb ions are less sensitive to strain^. Since PbTiC>3 
also has large atomic displacements it is likely to exhibit 
similar behaviour to PZT. 

The magnitudes of the calculated polarization are sim- 
ilar to those observed in the (100) epitaxial case for both 
perovskites 2 ^, apart from our predicted suppression of 
BaTiC>3 polarization at compressive strain. Such sup- 
pression under compressive strain has also been observed 
in (111) epitaxial BaTiO ^ 27 ' 28 , with remanent polariza- 
tion of only 0.1 ^Ccm~ 2 , opposed to 19.9 ^Ccm~ 2 in 
the (100) case^ 7 -. We are aware of no other reports on 
the magnitude of observed polarization in (111) epitaxial 
PbTi0 3 or BaTi0 3 . 

The calculations presented here do have a number 
of limitations. First, the calculations give the internal 
energy at zero temperature. A full strain-temperature 
phase diagram of a (111) epitaxial perovskite could be 
produced by parameterising an effective Hamiltonian 
model using first-principles results^. Second, the effect 
of the surface has not been considered. Third, no cell- 
doubling phases, such as octahedral rotations or other 
complex structures that require larger than 5-atom prim- 
itive cells, were considered, even though they may appear 
in perovskites in the stress regime considered here (e.g. 
Ref. [3]]) • The stability of the phases suggested here could 
be confirmed, for example, by calculating phonon modes. 
This was left beyond the scope of this study due to the 
large computational requirements of investigating each 
phase under various strains. Fourth, ionic zero-point mo- 
tion has a significant effect on the adopted phases as a 
function of pressure at low temperatures, and it could 
significantly alter the phase diagram^. For BaTiC>3 in 
particular, the shallowness of the potential wells makes 
this effect considerable up to several hundred Kelvin in 
the GPa pressure ranged 2 -, corresponding to strains con- 
sidered here. This effect could be taken into account 
by an effective Hamiltonian, if a path-integral quantum 
Monte Carlo technique, instead of classical Monte Carlo, 
were used in thermodynamic simulations. A final point to 



be noted is that while LDA is known to well reproduce 
many perovskite properties, the fact that it somewhat 
underestimates atomic displacements when compared to 
experiments could have an effect on the energetics of the 
various phases, as their elastic enthalpies are very close 
to each other. 
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